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Abstract 

An exact analytical solution of the statistical multifragmentation model is found in 
thermodynamic limit. The model exhibits a 1-st order phase transition of the liquid-gas 
type. The mixed phase region of the phase diagram, where the gas of nuclear fragments 
coexists with the infinite liquid condensate, is unambiguously identified. The peculiar 
thermodynamic properties of the model near the boundary between the mixed phase 
and the pure gaseous phase are studied. The results for the caloric curve and specific 
heat are presented and a physical picture of the nuclear liquid-gas phase transition is 
clarified. 
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Nuclear multifragmentation is one of the most interesting and widely discussed phenomena 
in intermediate energy nuclear reactions. The statistical multifragmentation model (SMM) 
(see H and references therein) was recently applied to study the liquid-gas phase transition 
in nuclear matter ^, |^. Numerical calculations within the canonical ensemble exhibited 
many intriguing peculiarities of the finite multifragment systems. However, the investigation 
of the system's behavior in the thermodynamic limit was still missing. Therefore, there was no 
rigorous proof of the phase transition existence, and the phase diagram structure of the SMM 
remained unclear. Previous numerical studies for the finite nuclear systems (the canonical and 
microcanonical ensembles) led to unjustified (and sometimes wrong) statements concerning the 
nuclear liquid-gas phase transition in the thermodynamic limit. In our recent paper ^ an 
exact analytical solution of the SMM was found within the grand canonical ensemble which 
naturally allowed to study the thermodynamic limit. The self-consistent treatment of the 
excluded volume effects was an important part of this study. In this letter we investigate the 
peculiar thermodynamic properties near the boundary between the mixed phase and the pure 
gaseous phase. New results for the caloric curve and the specific heat are presented and a 
physical picture of the nuclear liquid-gas phase transition in SMM is clarified. This physical 
picture differs from the one advocated in the previous numerical studies. 

In the SMM the states of the system are specified by the multiplicity sets {uk} {uk = 
0, 1, 2, ...) of fc-nucleon fragments. The partition function of a single fragment with k nucleons 
is [jl|: LUk = V {mT k / 271)^ ^"^ , where = 1, 2, A [A is the total number of nucleons in the 
system). V and T are, respectively, the volume and the temperature of the system, m is the 
nucleon mass. The first two factors in uo^ originate from the non-relativistic thermal motion 
and the last factor, Zk-, represents the intrinsic partition function of the /c-fragment. For k = 1 
(nucleon) we take zi = 4 {4 internal spin-isospin states) and for fragments with > 1 we use 
the expression motivated by the liquid drop model (see details in Ref. |^): z^ = exp(— /^/T), 
with the fragment free energy 

fk = - [Wo + TVeo ]k + a{T) + r Tlnfc . (1) 

Here Wo = 16 MeV is the bulk binding energy per nucleon, T^/eo is the contribution of 
the excited states taken in the Fermi-gas approximation (eo = 16 MeV) and cr(T) is the 
temperature dependent surface tension which is parameterized in the following form: o"(T) = 
cTo[(T2 - T2)/(T2 + T2)]5/4, with do = 18 MeV and T, = 18 MeV (a = at T > T,). 
The last Fisher's term in Eq. (0) with dimensionless parameter r is introduced for generality. 
The canonical partition function (CPF) of the ensemble of nuclear fragments has the following 
form: 

Z'^{V,T) = E n '^SiA-Y.knk) . (2) 

The model defined by Eqs.(|l],|2D with r = was studied numerically in Refs. H, ^, |[. This 
is a simplified version of the SMM since the symmetry-energy and Coulomb contributions 
are neglected. However, its investigation appears to be very important for understanding the 
physics of multifragmentation. 

In Eq. (^ the nuclear fragments are treated as point-like objects. However, these fragments 
have non-zero proper volumes and they should not overlap in the coordinate space. In the Van 
der Waals excluded volume approximation this is achieved by replacing the total volume V in 
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Eq. (0) by the free (available) volume Vf = V—bJ2k knk, where b = 1/po (Po = 0.16 fm~'^ is the 
normal nuclear density). Therefore, the corrected CPF becomes: Za{V,T) = Z'^^{V — bA,T). 

The calculation of Za^V-, T) is difficult because of the constraint J2k kuk = A. This difficulty 
can be partly avoided by calculating the grand canonical partition function: 

oo 

Z{V, T, /i) = exp (M/T) ZAiV, T) Q{V - bA) , (3) 

where the chemical potential /i is introduced. The calculation of Z is still rather difficult. 
The summation over the sets {uk} in Za cannot be performed analytically because of the 
additional A-dependence in the free volume Vf and the restriction V/ > 0. The problem can 
be solved by introducing the so-called isobaric partition function (IPF) which is calculated in 
a straightforward way (see details in Refs. [0, @]): 

Z(s,T,/x) ^ r dV eM-sV) Z{V,T,^) = ^) ^ . , (4) 

JO s — j^isU.u) 



where 



J'(s,T, /i) 

fi-sbT\ , ^^y^_^_f{u-sbT)k-a0^' 

k=2 



zi exp I '-^ ) + J2 ' 



, (5) 



with u = fi + Wo + T^/eo. In the thermodynamic limit V ^ oo the pressure of the system is 
defined by the farthest-right singularity, s*{T,fi), of the IPF Z{s,T,fi) 

p{T,^,) ^ T hm ^^^'^^ = T .*(r,/.) . (6) 

V -^oo y 

The study of the system's behavior in the thermodynamic limit is therefore reduced to the 
investigation of the singularities of Z. 

The IPF (^) has two types of singularities: 1) the simple pole singularity defined by 
the equation Sg(T, /i) = J^{sg,T,fi) ; 2) the singularity of the function itself at the point 
s;(T, /i) = u/Tb where the coefficient in linear over k terms of the exponent in Eq. (^ is equal 
to zero. 

The simple pole singularity corresponds to the gaseous phase where pressure pg = Tsg 
is determined by the transcendental equation: pg(T,fi) = TJ^{pg/T,T, ft). The singularity 
si(T,fi) of the function defines the liquid pressure: pi(T,^) = Tsi(T,fi) = u/b. Here the 
liquid is represented by an infinite fragment (condensate) with k = oo. 

In the region of the (T, /i)-plane where z/ < bpg{T, fi) the gaseous phase is realized {pg > pi), 
while the liquid phase dominates at > bpg{T, /i). The liquid-gas phase transition occurs when 
the two singularities coincide, i.e. Sg{T,fi) = si{T,fi). As JF in Eq. (|) is a monotonously 
decreasing function of s the necessary condition for the phase transition is that the function 
J-' is finite in its singular point si. At r = this condition requires ct(T) > and, therefore, 
T < Tc. Otherwise, JF(s;, T, fi) = oo and the system is always in the gaseous phase as Sg > si. 
As one can see from Eq.(^ the convergence properties of J-'{s,T,^) depend significantly on 
the Fisher's exponent r in the vicinity of the critical point where the surface term vanishes. In 
what follows we concentrate on the case r = 0. Other possibilities which appear at r > are 
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discussed in Ref. |^ and their detail study will be presented elsewhere. Here we only note that 
Eqs. ^ represent an exact solution of the Fisher's droplet model where additionally 
the effects of excluded volume are incorporated. 

The baryonic density p in the liquid and gaseous phases is given by the following formulae, 
respectively: 

Pi = idpi/dfi)j, = 1/6 , pg = {dpg/dp)j, = pid/{l + hpid) , 

where the function pid is the density of point-like nuclear fragments with shifted, p ^ p — bpg, 
chemical potential: 



p.d{T,p)=[^—) 



(7) 



A similar expression for pg within the excluded volume model for the pure nucleon gas was 



obtained in Ref. 10 



At T < Tc the system undergoes a 1-st order phase transition across the line p* = p*{T) 
defined by the condition of coinciding singularities: si = Sg, i.e., pi = Pg. The phase transition 
line p*(T) in the (T, /i)-plane corresponds to the mixed liquid and gas states. This line 
is transformed into the finite mixed-phase region in the (T, p)-plane shown in Fig. 1. The 
baryonic density in the mixed phase is a superposition of the liquid and gas baryonic densities: 
p = Xpi + {1 — X)pg , where A (0 < A < 1) is the fraction of the system's volume occupied by the 
liquid inside the mixed phase. Similar linear combinations are also valid for the entropy density 
s and the energy density e with {i = l,g) Si = (dpi/dT)^ , Si = T (dpi/dT) ^ + p {dpi/dp)rp — pi. 

Inside the mixed phase at constant density p the parameter A has a specific temperature 
dependence shown in Fig. 2: from an approximately constant value p/po at small T the 
function A(T) drops to zero in a narrow vicinity of the boundary separating the mixed phase 
and the pure gaseous phase. This corresponds to a fast change of the configurations from the 
state which is dominated by one infinite liquid fragment to the gaseous multifragment config- 
urations. It happens inside the mixed phase without discontinuities in the thermodynamical 
functions. 

An abrupt decrease of A(T) near this boundary causes a strong increase of the energy 
density as a function of temperature. This is evident from Fig. 3 which shows the caloric curves 
at different baryonic densities. One can clearly see a leveling of temperature at energies per 
nucleon between 10 and 20 MeV. As a consequence this leads to a sharp peak in the specific 
heat per nucleon at constant density, Cp(T) = {de/dT)^/ p , presented in Fig. 4. A finite 
discontinuity of Cp(T) arises at the boundary between the mixed phase and the gaseous phase. 
This finite discontinuity is caused by the fact that A(T) = 0, but {d\/dT)p 7^ at this 
boundary (see Fig. 2). 

It should be emphasized that the energy density is continuous at the boundary of the mixed 
phase and the gaseous phase, hence the sharpness of the peak in Cp is entirely due to the strong 
temperature dependence of A(T) near this boundary. Moreover, at any p < po the maximum 
value of Cp remains finite and the peak width in Cp(T) is nonzero in the thermodynamic limit 
considered in our study. This is in contradiction with the expectation of Refs. Q that 
an infinite peak of zero width will appear in Cp(T) in this limit. Also a comment about the 
so-called "boiling point" is appropriate here. This is a discontinuity in the energy density 
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£{T) or, equivalently, a plateau in the temperature as a function of the excitation energy. Our 
analysis shows that this type of behavior indeed happens at constant pressure, but not at 
constant density! This is similar to the usual picture of a liquid-gas phase transition. In Refs. 
0, ^ a rapid increase of the energy density as a function of temperature at fixed p near the 
boundary of the mixed and gaseous phases (see Fig. 3) was misinterpreted as a manifestation 
of the "boiling point" . 

In conclusion, the simplified version of the SMM is solved analytically in the grand canon- 
ical ensemble. The progress is achieved by reducing the description of phase transitions to the 
investigation of the isobaric partition function singularities. The model clearly demonstrates 
a 1-st order phase transition of the liquid-gas type. The considered system has peculiar 
properties near the boundary of the mixed and gaseous phases. The rapid change of the 
thermodynamical functions with T at fixed p takes place near this boundary due to the dis- 
appearance of the infinite liquid fragment. This leads to leveling of the caloric curves shown 
in Fig. 3 at temperatures between 6-10 MeV depending on the density. As a consequence a 
sharp peak and a finite discontinuity are developed in the specific heat Cp{T) at the boundary 
of the mixed and gaseous phases. 
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Figure 1: Phase diagram in the (T, p)-plane. The mixed phase and pure gaseous phase 
boundary is shown by the sohd hne. The pure hquid phase (shown by crosses) corresponds to 
the fixed density p — po- Point C is the critical point, at T > Tc only the pure gaseous phase 
exists. 




T(MeV) 

Figure 2: Volume fraction A(T) of the liquid inside the mixed phase is shown as a function 
of temperature for fixed nucleon densities p/po — 1/6, 1/3, 1/2, 2/3, 5/6 (from bottom to top). 
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Figure 3: Temperature as a function of energy density per nucleon (caloric curve) is shown 
for fixed nucleon densities p/po — 1/6, 1/3, 1/2, 2/3. 
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Figure 4: Specific heat per nucleon as a function of temperature at fixed nucleon density 
p/po — 1/3. The dashed line shows the finite discontinuity of Cp{T) at the boundary of the 
mixed and gaseous phases. 
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